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Abstract 

in ■ 

We have searched for absorption in several common atomic species due to the atmosphere or exospherc 
of the transiting extrasolar planet HD 209458b, using high precision optical spectra obtained with the 
Subaru High Dispersion Spectrograph (HDS). Previously we reported an upper limit on Ha absorption of 
0.1% (3c) within a 5.1 A band. Using the same procedure, we now report upper limits on absorption due 
■ to the optical transitions of Na D, Li, Ha, H/3, H7, Fe, and Ca. The 3cr upper limit for each transition is 

Oh I approximately 1% within a 0.3 A band (the core of the line), and a few tenths of a per cent within a 2 A 

q ■ band (the full line width). The wide-band results are close to the expected limit due to photon-counting 

I (Poisson) statistics, although in the narrow-band case we have encountered unexplained systematic errors 

at a few times the Poisson level. These results are consistent with all previously reported detections and 
upper limits, but are significantly more sensitive. 
, Key words: stars: planetary systems: individual (HD 209458b) — techniques: spectroscopic 

?-h I 1. Introduction the starlight is transmitted through the thin, partially- 
transparent portion of the planetary atmosphere and ex- 

The first star other than the Sun for which a planetary osphcre. Thus, by obtaining spectra with high spectral 

transit was detected was HD 209458. The initial discov- resolution and high signal-to-noise ratio, and comparing 

ery by Charbonneau et al. (2000) and Henry et al. (2000), the spectra taken in and out of transit, one can detect 

along with subsequent observations of higher accuracy by absorption or emission features due to the planetary at- 

Brown et al. (2001) and Jha et al. (2000), allowed mea- mosphere or exosphere. This methodology is referred to 

suremcnts of the orbital inclination (i = 86?1) and the as "transmission spectroscopy." 

planetary radius (R — 1.34 i?j up ). By combining these Several studies of HD 209458 using transmission spec- 
results with the stellar radial velocity data, it was shown troscopy have previously been reported. Most have pro- 
that the companion has a mass comparable to Jupiter vided upper limits on particular absorption features. The 
(M = 0.69 Mj U p) and the density of a gas giant planet first positive detection of an extrasolar planetary atmo- 
(p ~ 0.4 g cm -3 ). spheric absorption feature was made by Charbonneau 
Observations of planetary transits also provide an op- et al. (2002). Their motivation was to confirm the 
portunity to study the atmospheric properties of close-in prediction by Seager and Sasselov (2000) (and later, 
giant extrasolar planets (also known as "hot Jupiters"). Brown 2001; Hubbard et al. 2001) of a relatively strong 
During occultations, a small fraction of the stellar disk absorption feature in the sodium doublet lines. They 
is blocked from view, and an even smaller fraction of observed HD 209458 with the Space Telescope Imaging 
* ~ Spectrograph (STIS) onboard the Hubble Space Telescope 
U c c ow - (HST), and found an additional 0.023 ±0.006% absorp- 
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tion during transits within a band pass centered on the 
sodium doublet lines (5893 ± 6 A). However, this addi- 
tional absorption was significantly weaker than that pre- 
dicted. Later, detection of excess absorption in the Lya 
(1216 A) transition of neutral hydrogen was reported by 
Vidal-Madjar et al. (2003), also based on HST/STIS spec- 
tra. These authors found an amazingly strong additional 
Lya absorption of 15 ±4% during transits, implying that 
the absorbing H I gas extends over several planetary radii, 
and thus should be interpreted as an exosphere rather 
than an atmosphere. Vidal-Madjar et al. (2004) subse- 
quently reported 2-3 a detections of oxygen (O i) and 
carbon (C II) using HST/STIS. Interestingly, among the 
seven known cases of transiting hot Jupiters (Udalski 
et al. 2002a; Udalski et al. 2002b; Udalski et al. 2003; 
Alonso et al. 2004), HD 209458b has both the smallest 
mass and the largest radius, and thus an anomalously low 
density. This fact makes it yet more interesting to study 
the very extended exosphere of HD 209458b. 

These reports of an extended exosphere around 
HD 209458b have inspired searches for additional absorp- 
tion features. However, in spite of the successful detec- 
tions via space-based observations, no data obtained from 
the ground has yet confirmed the detections or discov- 
ered any new components. Upper limits on planetary ex- 
osphcric absorption features of various components in the 
optical region were reported by Bundy and Marcy (2000) 
based on Keck I/HIRES data. They obtained transmis- 
sion spectra of HD 209458 (H/3, H 7 , Fc) and 51 Pegasi 
(Li, Ha, Na Dl, and D2), and placed upper limits of 3- 
20% on additional absorption during transits in a 0.3 A 
bandwidth centered on the cores of the relevant stellar 
lines. Moutou et al. (2001) also searched for transit ef- 
fects (especially for ionized species) in the optical band 
with the VLT/UVES. They did not set upper limits on 
any specific atomic lines, but estimated a detection limit 
of order 1% within any 0.2 A band covered by the observa- 
tions. Near- infrared features have also been examined to 
set more stringent upper limits on He, CO, H2O, and CH4 
(e.g. Brown et al. 2002; Harrington et al. 2002; Moutou 
et al. 2003; Doming et al. 2005). 

In this paper we report on the results of searches for 
excess Na (Dl, D2), Li, H/3, Hy, Fe and Ca absorption 
in optical spectra of HD 209458 during transits, based 
on Subaru High Dispersion Spectrograph (HDS) data ac- 
quired in 2002 October. Compared with the previous in- 
vestigations, the HDS data cover a larger range of wave- 
lengths (covering the entire optical band) and have a 
higher signal-to-noise ratio. Moreover, this data set has a 
substantial advantage of covering a wide range of orbital 
phases, namely before, during and after the transit in a 
single night. The major obstacles to ground-based trans- 
mission spectroscopy are time-dependent effects due either 
to instrumental changes caused by variations in position, 
temperature et cetera, or to telluric effects such as seeing, 
airmass, and atmospheric absorption. All such effects are 
most effectively monitored, interpolated and removed if 
the data are obtained on a single night and if an entire 
transit is observed. We employed the same reduction and 



analysis techniques used and described in detail in Winn 
et al. (2004) (hereafter Paper I) in which we set an upper 
limit of 0.1% (3cr) additional absorption in a 5.1 A band 
centered on Ha. 

This paper is organized as follows. Section 2 describes 
the observations, summarizes the reduction method and 
points out some small differences between the techniques 
employed in the present paper and Paper I. We present 
results for the selected absorption lines in section 3 and 
discuss systematic effects in section 4. Section 5 provides 
a brief summary and interpretation of the results. 

2. Observations and Reduction Method 

We observed HD 209458 on the nights of 
HST 2002 October 24 and 26 (25 and 27 in UT) 
with the Subaru HDS (Noguchi et al. 2002). The planet 
was predicted to transit its parent star near the middle 
of the night of HST October 24 (see figure 1), according 
to the ephemcris of Brown et al. (2001). In order to 
determine accurate radial velocities of the star and to 
verify the ephemeris, we obtained a spectrum with an 
Iodine Cell at both the start and the end of this series of 
exposures, and confirmed our expectation that the transit 
took place when predicted. We obtained 30 spectra of 
HD 209458, and 7 spectra of the rapidly rotating B5 Vn 
star HD 42545 (Abt et al. 2002) in order to evaluate the 
telluric spectrum. The wavelength coverage was 4100 A 
< A < 6800 A with a spectral resolution of R m 45000 
(0.9 km s _1 per pixel). The typical exposure time was 
500 seconds with seeing of ~ 0'.'6, through air masses 
ranging from 1.0 to 1.9. The resulting spectra had a 
typical signal-to-noise ratio of ~ 350 per pixel. 
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First, we processed the frames with standard IRAF 1 
procedures and extracted one-dimensional spectra. Next, 
we applied an analysis method described at length in 
Paper I to correct for substantial time-dependent varia- 
tions of the instrumental blaze function. However, a slight 
modification of the Paper I procedure was employed. In 
Paper I, we used all 4100 pixels of each spectral order to 
calculate the best-fitting parameters of a smooth function 
that describes the instrumental variations; however, we 
noted that this procedure acts to dilute any real signal. 
An artificially injected signal of 0.1% absorption was re- 
covered with a strength of only 0.057%. We have since 
found that this problem is significantly reduced if the do- 
main of pixels used in the fitting procedure is restricted 
to approximately 1000 pixels in the vicinity of the tar- 
get line. This modification increases the sensitivity of the 
analysis and also provides a modestly improved correction 
for instrumental variations (see section 4). 

Next, a template spectrum T(A) is created from all 30 
individual exposures using the above correction method 
(see Paper I for further details). A time series of residual 
spectra are generated by subtracting the suitably matched 
template T m (A) from each spectrum S(X,t); thus the tem- 
plate spectrum determines the zero-level for the residu- 
als, and therefore we use all of the 30 spectra (including 
in-transit spectra) so as to create an optimum template 
spectrum. Then, we examine the time variation of each 
selected absorption line in order to search for possible ad- 
ditional absorption in the transit phase. Specifically, we 
compute a time series of fractional difference in flux be- 
tween the spectrum and the matched template within a 
smoothing width AA according to 



5 A (t) 



Ao+AA/2 



Ao-AA/2 



dX [S(X,t)-T m (X)} 



(1) 



™A +AA/2 ' 

dX T m (X) 

'Ao-AA/2 

where A indicates each target absorption line. We refer 
to such a time variation of the fractional difference as a 
"difference light curve" . 

To place a quantitative limit on additional absorption, 
we calculate <5; n a and a- m a- 
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which are the mean value and the standard deviation of 
i5a (i) for the spectra taken in the transit (between sec- 
ond and third contacts, namely within ± 60 minutes from 
center of the transit). Similarly, we compute <J ut,A and 
Cout,A from the spectra taken out of the transit (before 

The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the U.S. National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



first contact or after fourth contact: beyond ± 92 minutes 
from center of the transit). 2 The numbers of spectra in 
each phase are N- in = AT out = 12. 
Finally, we compute 



5 a = Sint 



Jout.A i 



Noutcr 2 



CA 



out, A 



'out A 



N in + N a 



(4) 



(5) 



Note that 5a would be negative if there were additional 
absorption during transit. 

In the following sections, we express the results as 5a ± 
a a- Where no additional absorption feature related to 
transit are detected (i.e., \5a \ < o~a), we quote upper limits 
as 3(Ta- 

3. Results 

Since the HDS spectra cover a wide wavelength region in 
optical band (4100 6800 A), there are many absorp- 
tion lines available for inspection. We concentrate on the 
following transitions for detailed study: 1) the sodium and 
lithium lines, because additional absorptions at these lines 
are theoretically predicted by Seager and Sassclov (2000); 
2) neutral hydrogen lines including Ha, H/3, and H7, since 
an extended hydrogen exosphere was reported by Vidal- 
Madjar et al. (2003); and 3) the strong stellar lines of 
Ca and Fe. Another reason to examine the sodium dou- 
blet lines in detail is to attempt to confirm the report by 
Charbonneau et al. (2002) of additional 0.023 ± 0.006% 
absorption in the 5893 ± 6 A band. 

In Paper I, we paid particular attention to the results for 
a bandwidth of 5.1 A, in order to facilitate a comparison 
with the band used by Vidal-Madjar et al. (2003). In this 
work we generally consider two different bandwidths: a 
relatively narrow band of 0.3 A [so as to facilitate compar- 
ison with the previous results of Bundy and Marcy (2000) 
and Moutou et al. (2001)], and a relatively wide band of 
2 A to cover the full width of target absorption lines and 
provide more sensitivity. 

3.1. Calculation of Difference Spectra 

We express the difference between each spectrum and 
the matched template in units of the expected standard 
deviation due to Poisson statistics, 



S(X)-T m (X) 

\/?m(A) 



(6) 



and refer to this quantity as a "difference spectrum" . As 
an example, the difference spectra in the vicinity of the 
sodium resonance doublet are shown in figure 2. The tem- 
plate spectrum and the telluric spectrum are plotted at 
the bottom of the series. To estimate the telluric spec- 
trum, we created a template spectrum for HD 42545 in 



2 To compute the transit times, we use an ephemeris defined by 
J.D. (mid-transit) = 2451659.93675 (Brown et al. 2001) and the 
most precise published period, P = 3.524739 d (Robichon and 
Arenou 2000). 
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5880 5885 5890 5895 5900 

wavelength [Angstroms] 

Fig. 2. Difference between each spectrum and its matched template for a wavelength range in the vicinity of the sodium doublet 
lines. The spectra are ordered in time from bottom to top. The corresponding phases of the planetary transit are noted on the 
right side of the figure. The spacing between minor ticks represents two Poisson deviations. The bottom spectrum shows the 
normalized template spectrum of the rapidly rotating star HD 42545, containing telluric features and interstellar sodium absorption. 
The second spectrum shows the overall template spectrum of HD 209458. Note that we do not correct to the laboratory frame. 



the same way as for HD 209458 and normalized it appro- 
priately for visual reference. 

The template spectrum shows some features that are 
similar to those seen in the estimated telluric spec- 
trum. However, the latter spectrum is in reality not 
just the telluric spectrum; it is also affected by interstel- 
lar sodium absorption. Moreover, the sodium D2 line at 
5890 A is blended with strong telluric vapor absorption 
(e.g., Lundstrom et al. 1991), and the variation of meso- 
spheric sodium column density is not negligible (e.g., Ge 
et al. 1997). These factors substantially complicate any 
interpretation of the observed Na lines. 

There are time-dependent, coherent features with am- 
plitudes that are in excess of the Poisson noise level in 
the vicinity of the Na lines (such features are also seen for 



other lines, but at a much reduced level). These anoma- 
lies strongly resemble those which would be produced by 
a small, time-dependent wavelength shift. However, mea- 
surements of the line center wavelengths (obtained by fit- 
ting a Gaussian function to the line profile) of both stellar 
and telluric features show them to be stable to within 
±0.01 A (less than 1 pixel), and show no systematic drift 
with time. Thus, they do not appear to be directly asso- 
ciated with errors or drift in our wavelength calibration 
(determined from Thorium-Argon lamp spectra obtained 
at the beginning and end of the HD 209458 exposures). 
Furthermore, these features are strongest outside of the 
transit phase, especially at the beginning of the observa- 
tion when the target was near zenith, and thus when the 
instrument was rotating most rapidly. Hence, we believe 



No. 



Transmission Spectroscopy of the HD 209458b with Subaru HDS 



ii 

< 



I 

CO 



0.01 



-0.01 



6 0.01 



-0.01 







-0.01 7 



0.01 - 







~i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

D1 (blue side) 



i 
i 



5895.06±0.15A 

I I I I I I I I I I 



-I— H 



Mi' 
1 1 i 



D1 (line core) 



-I 



ill i 



i 



I Si I if 



i 



5896.06±0.15A 



I I I I | I I I I 
D1 (red side) 



-I— H 



xL 

III 1 2 2 f 
I 
I 



I 



1*1 



5897.06±0.15A 

J i i i i I i i i i L 



n — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

D2 (blue side) 



i 
i 



5889.09±0.15A 



I NI I I 



I I I I 



-Hi 



I 



D2 (line core) 



Ill ^1 



I 



5890.09±0.15A 



I NI I I 



I I I I 



-I— H 



D2 (red side) 

I ; I 

I | I 

I : : I 

I : : I 

5891.09±0.15A 

i i i i I i i i i I i i i i I i 



■100 







100 

time [min] 



-100 







100 



Fig. 3. Difference light curves around the sodium doublet lines for a smoothing width of AA = 0.3 A. The error bars repre- 
sent the Poisson noise \/7m, and the horizontal axis indicates the time of each exposure with the origin at mid-transit. The 
time between the two dotted lines indicates the in-transit phase, and those between the dashed line and dotted line show 
ingress (left side) and egress (right side) phases. The left side panels exhibit difference light curves for the Na Dl line at 
5895.06 ±0.15 A (top panel; continuum of blue side), 5896.06 ±0.15 A (middle panel; line core), 5897.06 ±0.15 A (bottom 
panel; continuum of red side). The right side panels are similar, but for the Na D2 line; 5889.09 ±0.15 A (top panel; con- 
tinuum of blue side), 5890.09 ±0.15 A (middle panel; line core), 5891.09 ±0.15 A (bottom panel; continuum of red side). 



the anomalies are unlikely to be due to real absorption in 
the HD 209458 system. Whatever their source, these fea- 
tures appear to have almost exactly zero net photon count; 
they cancel out almost completely when integrated over 
a smoothing width comparable to the overall line width 
(e.g., 2 A). This suggests that they do not compromise 
the achieved detection sensitivity for such a bandwidth. 

3.2. Comparison with Previous Results 

The same analysis was applied to the other selected 
lines. For the other lines, there was no significant in- 
terstellar absorption contamination such as that observed 
for Na. In order to compare our results with previous 



ones, we selected a bandwidth matched to past studies. 
In particular we followed the conventions of Bundy and 
Marcy (2000), who performed transmission spectroscopy 
for HD 209458 and 51 Pegasi in various atomic lines. They 
used AA ~ 0.3 A in order to emphasize the line cores. 

The time variations of the fractional difference for AA = 
0.3 A at the sodium doublet lines are plotted in figure 3. 
It is evident that the difference light curves for contin- 
uum regions are nearly consistent with the Poisson noise 
of ~ 0.07%, while those for line cores have a significantly 
larger scatter than the expected Poisson noise of ~ 0.15%. 
The results for line cores are <5di = 0.12 ± 0.48% and 
<!>D2 = —0.70 ±0.28%, respectively. Although (5d2 indicates 
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Fig. 4. Comparison of the Subaru HDS upper limits for AA = 0.3 A with those of Bundy and Marcy (2000). The ver- 
tical axis is fractional upper limits on additional absorption during transit on a log scale. The horizontal axis indi- 
cates atomic species and corresponding line wavelength Ao in equation (1). Note that Bundy and Marcy (2000) did 
not report a Ca (4226.82 A) result, and their upper limits on Li, Hce, and Na (Dl, D2) are for 51 Pegasi. 
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Fig. 5. Difference light curves for AA = 2 A around the sodium doublet. The left panels exhibit the 

light curves on the transiting night (the night of HST 2002 October 24), and the right panels are 
those for the out-of-transit night (October 26). The vertical scale is matched to that of figure 3. 



No. ] 

a 2.5a detection of additional absorption, <5di does not ex- 
hibit a corresponding signal, and as noted at the end of 
the previous subsection, the difference spectra have unex- 
plained systematic anomalies on a wavelength scale larger 
than the 0.3 A smoothing region. Thus, we do not regard 
this 2.5a signal as real absorption, but rather as an indi- 
cator of systematic error. The resulting upper limits are 
3ct D i = 1-43% and 3cr D2 = 0.85%. 

Figure 4 shows a comparison of our upper limits for 
HD 209458 with those of Bundy and Marcy (2000) for 
HD 209458 (H/3, H 7 , Fe) and 51 Pcgasi (Li, Ha, Na Dl, 
and D2). Note that although 51 Pcgasi does not show 
transits, they consider that the lack of transits does not 
rule out the possibility of excess absorption features. The 
Subaru HDS upper bounds are the most stringent limits 
so far that have been obtained with ground-based obser- 
vations in the optical region. 

3.3. Results for Integration over the Full Line Width 

We repeated the same procedure described in the pre- 
vious subsection using a smoothing width A A = 2 A in 
order to cover the typical width of the strong stellar ab- 
sorption lines. In this case, the contribution from the 
small systematic features discussed in subsection 3.1 was 
nearly zero, because of the inversion symmetry of those 
anomalies. Wc display the difference light curves of the 
Na D lines in figure 5, as examples. The large scatter 
seen in figure 3 is almost entirely removed for this larger 
smoothing width. The results are foi = 0.06 ± 0.09% and 
8-di = -0.02 ± 0.06%. We also plot S(t) for the out-of- 
transit data obtained on the night of HST 2002 October 
26 for reference (right side). These panels exhibit very 
similar variations to those in the October 24 data, and 
thus support the conclusion of no detected additional ab- 
sorption during the transit event. 

Our fractional upper limits on additional absorption 
during transit for all of the selected lines and both val- 
ues of A A lines are given in table 1. We did not set upper 
limits on Fe I and Ca I in blue CCD, since the bands of 
AA = 2 A contain several irrelevant absorption lines. 

We also derived an upper limit on the Na D lines 
with AA = 12 A, corresponding to the band in which 
Charbonneau et al. (2002) reported detection of addi- 
tional absorption (5893 ±6 A). The resulting difference 
light curve is plotted in figure 6, and the result is 5-^^ = 
0.03 ±0.04%. Thus, we are not able to either confirm or 
contradict the additional absorption of 0.02% reported by 
Charbonneau et al. (2002), since it is beyond our sensitiv- 
ity level. 

4. Evaluation of Systematic Errors 

In order to evaluate the significance of high precision 
measurements, such as those presented in this paper and 
to facilitate the planning of similar, future observations, 
it is necessary to thoroughly understand the sources and 
magnitudes of systematic errors. 

Wc therefore carried out further analyses of the remain- 
ing systematic effects, such as instrumental variations and 
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Table 1. Summary of the 3cr upper limits. 



Element 


Ao 
(A) 


AA = 0.3 A 

(%) 


AA = 2 A 

(%) 


red CCD 


Li I 
Ha 
Fe I 

Na i(Dl) 
Na i(D2) 


6707.95 
6562.95 
6024.20 
5896.06 
5890.09 


0.91 
0.81 
0.84 
1.43 
0.85 


0.43 
0.66 
0.18 
0.28 
0.18 


blue CCD 


TT O 

tip 
Fe I 

H 7 
Ca i 


4861.45 
4383.65 
4340.56 
4226.82 


0.68 
1.08 
1.57 
2.14 


0.34 
0.88 
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Fig. 6. Difference light curve for 5893 ±6 A, correspond- 
ing to the band in which Charbonneau et al. (2002) re- 
ported the detection of 0.02% additional absorption. 

telluric effects, which may be serious obstacles to improv- 
ing ground-based observations. First of all, we checked the 
reliability of our correction method for the instrumental 
variations and estimated the magnitude of any remaining 
instrumental effects on the results. Next, we confirmed 
that the correction method does not dilute any real signal 
which is above the Poisson noise limit. Finally, we ana- 
lyzed the correlation between the spectra of HD 209458 
and rapidly rotating star HD 42545 in order to estimate 
the influence of terrestrial atmospheric effects. 

4-.1. Estimation of Instrumental Effects 

Although our spectrum-matching procedure (Paper I 
gives its details) is intended to correct for instrumental 
variations, the correction is not perfect. To estimate the 
size of the remaining instrumental effects, we selected an 
appropriate reference absorption line on the red CCD. We 
adopted the Fc i] (intercombination) line at 6024.20 A for 
the following reasons: (i) there are no telluric lines within 
±2 A (Moore et al. 1966), and (ii) the position of the line 
is near the peak of the continuum, which minimizes the 
Poisson noise level. The same procedure described in the 
previous section was applied with AA = 0.3 A and 2 A to 
the band around 6024.20 A. 
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Fig. 7. Difference light curves for the reference Fe line at 6024 A. The upper panels show frac- 
tional differences of AA = 0.3 A, while the lower panels use AA = 2 A. The standard deviations 
are 0.28% (upper left), 0.65% (upper right), 0.06% (lower left), and 0.10% (lower right) respectively. 
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Fig. 8. Difference light curves with artificial absorption (left panels). The five data points between the dot-dashed lines in- 
dicate the results with artificial signals. The left panels show its recovery via our analysis procedure, and the right pan- 
els show residuals that subtracted the original difference light curve (lower left panel of figure 7.) from the left panels. 
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The difference light curves are shown in figure 7, and 
have standard deviations of 0.28% (upper left), 0.65% (up- 
per right), 0.06% (lower left), and 0.10% (lower right). As 
suggested by the analysis presented in the previous sec- 
tions, instrumental variations appear to be removed ap- 
proximately within the Poisson noise limit for the larger 
value of AA that integrates over the whole line width. 
However, for the narrower band, which includes only the 
line core, systematic effects of unknown origin increase the 
effective noise to a level that is several times larger than 
the expected Poisson noise. 

4-2. Recovery of Artificial Input Signals 

In order to verify that our reduction and analysis pro- 
cedures do not remove or dilute real signals, we injected 
an artificial signal of 0.03% (corresponding to the Poisson 
noise) and 0.2% into 5 out of 30 spectra in the October 

24 data around 6024 ± 1 A, and analyzed it via the 
same methods. Figure 8 shows difference light curves 
for AA = 2 A (left side) and residuals that subtracted 
the original difference light curve from left panels (right 
side). We compared the mean and standard deviation of 

25 realizations without artificial injection (8 no ,o~ no ), and 5 
with artificial signals injected (<5art,Cart)- We find that the 
residual (<5 no — <5 art ) ± (v/25cr2 o /30 + 5cr2 rt /30) values are 
0.0296±0.0003% (originally 0.03%) and 0.1987±0.0009% 
(originally 0.2%). We thus conclude that the procedures 
used in this paper can recover small input signals over the 
Poisson noise limits with AA = 2 A quite satisfactorily and 
without significant dilution. 

4-3. Telluric Effects 

Spectra obtained from the ground naturally contain 
many telluric absorption features. We obtained spectra 
of the rapidly rotating star HD 42545 at the end of the 
night (airmass ~ 1.9) and used this data to identify many 
such lines listed in Moore et al. (1966). Many telluric wa- 
ter vapor and oxygen absorption lines are visible in the 
spectra of the red CCD, particularly. 

The difference light curves for one strong telluric oxy- 
gen line at 6287.75 A with both AA = 0.3 A (top panel) 
and AA = 2 A (middle panel) are shown in figure 9. 
The bottom panel shows the difference light curve for 
5885.98 ±0.15 A, which contains one of the strongest tel- 
luric water vapor lines near the Na D doublet, although 
it is 5 to 10 times weaker than the 6287.75 A oxygen line. 
These plots show that the telluric features are measured 
to roughly within the Poisson noise (as expected), and 
do not exhibit any anomalous behavior during the tran- 
sit phase, and vary most strongly at the end of the night 
when the air mass became considerable (~ 2.0). 

Based on these results we estimate the telluric water 
vapor contamination levels by comparison with the equiv- 
alent widths listed in Moore et al. (1966), and conclude 
that they were no more than 0.5% for AA = 0.3 A, and 
less than 0.1% (3<r) for AA = 2 A for all of the lines listed 
in table 1. 
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Fig. 9. Difference light curves of a telluric oxygen line 
(upper and middle) and a water vapor feature (lower; 
near the Na D doublet). The strengths of both 
features increased near the end of the night when 
the air mass became significantly larger with time. 



5. Summary and Conclusions 

An investigation of the optical spectra taken during a 
transit of the HD 209458 system with Subaru HDS re- 
vealed no additional absorption features associated with 
the transit. The typical achieved sensitivity level is ~ 1% 
(3cr) for AA = 0.3 A (the line cores) in fractional difference 
spectra, and a few times smaller for AA = 2 A(thc full line 
width). Although we cannot confirm the smaller level of 
absorption in the Na D line that was observed from space 
(Charbonneau et al. 2002), our upper limits for the other 
optical transitions are consistent with, but more stringent 
than, those previously achieved from the ground. 

Given the systematic effects that we encountered when 
working at very high precision, one may wonder whether 
ground-based observations can ever be competitive with 
space-based observations for this purpose. The answer is 
unclear, but there arc at least a few lines of future research 
that may help to improve the prospects for ground-based 
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observations. As we have argued, the primary obstacle 
to higher-precision differential spectroscopy seems to be 
small temporal variations in instrumental response; tel- 
luric contamination is not negligible, but is removable to 
some extent. Instruments that are purpose-built for in- 
creased stability (despite telescope motion and instrument 
rotation) may be required. Alternatively, it may be pos- 
sible to use the multi-object capability of some spectro- 
graphs to record simultaneous high-resolution, high-SNR 
spectra of additional objects besides the exoplanetary sys- 
tem, and perform differential spectrophotometry. Finally, 
an alternative method for probing exoplanetary atmo- 
spheres was recently proposed by Snellen (2004). The 
idea is to measure the line-dependence of the Rossiter- 
McLaughlin effect (e.g., Rossiter 1924; McLaughlin 1924; 
Queloz et al. 2000; Ohta et al. 2005), which is the spec- 
tral line distortion caused by the planet since it blocks 
a particular part of the rotation field of the face of the 
star. This method involves a comparison of the shapes of 
different spectral lines, rather than their depths relative 
to the continuum. Thus using this method in ground- 
based observations is quite reasonable, because ground- 
based instruments would not exceed space-based ones in 
photometric accuracy, but possibly could do so in spectral 
resolution. We intend to use the existing Subaru/HDS 
data to investigate this latter possibility. 
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